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ABSTRACT 

The MAGIC very-high-energy (VHE) 7-ray astronomy collaboration recently 
reported the detection of the quasar 3C279 at > 100 GeV 7-ray energies. Here we 
present simultaneous optical (BVRI) and X-ray [RXTE PCA) data from the day 
of the VHE detection and discuss the implications of the snap-shot spectral en- 
ergy distribution for jet models of blazars. A one- zone synchrotron-self-Compton 
origin of the entire SED, including the VHE 7-ray emission is highly problematic 
as it would require an unrealistically low magnetic field. The measured level 
of VHE emission could, in principle, be interpreted as Compton upscattering of 
external radiation (e.g., from the broad-line regions). However, such an interpre- 
tation would require either an unusually low magnetic field oi B ^ 0.03 G, or (in 
order to achieve approximate equipartition between magnetic field at i? ~ 0.25 G 
and relativistic electrons) an unrealistically high Doppler factor of P ~ 140. In 
addition, such a model fails to reproduce the observed X-ray flux. Furthermore, 
both versions of leptonic one-zone models produce intrinsic VHE 7-ray spectra 
steeper than measured, even in the case of the lowest plausible extragalactic 77 
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absorption. We therefore conclude that a simple one-zone, homogeneous lep- 
tonic jet model is not able to plausibly reproduce the SED of 3C279 including 
the recently detected VHE 7-ray emission. This as well as the lack of correlated 
variability in the optical with the VHE 7-ray emission and the substantial 77 
opacity of the BLR radiation field to VHE 7-rays suggests a multi-zone model 
in which the optical emission is produced in a different region than the VHE 
7-ray emission. In particular, an SSC model with an emission region far outside 
the BLR reproduces the simultaneous X-ray — VHE 7-ray spectrum of 3C279. 
Alternatively, a hadronic model is capable of reproducing the observed SED of 
3C279 reasonably well, both in scenarios in which only the internal synchrotron 
field serves as targets for jry pion production, and with a substantial contribu- 
tion from external photons, e.g., from the BLR. However, either version of the 
hadronic model requires a rather extreme jet power of up to Lj ~ 10^^ erg s~^, 
compared to a requirement of ~ 2 x lO"^^ erg s~^ for a multi-zone leptonic 
model. 

Subject headings: galaxies: active — Quasars: individual (3C 279) — gamma- 
rays: theory — radiation mechanisms: non-thermal 



1. Introduction 



Flat-spectrum radio quasars (FSRQs) and BL Lac objects are active galactic nuclei 
(AGNs) commonly unified in the class of blazars. They exhibit some of the most violent high- 
energy phenomena observed in AGNs to date. Their spectral energy distributions (SEDs) 
are characterized by non-thermal continuum spectra with a broad low-frequency component 
in the radio - UV or X-ray frequency range and a high-frequency component from X-rays to 
7-rays. In the framework of relativistic jet models, the low-frequency (radio - optical/UV) 
emission from blazars is interpreted as synchrotron emission from nonthermal electrons in a 
relativistic jet. The high-frequency (X-ray - 7-ray) emission could either be produced via 
Compton upscattering of low frequency radiation by the same electron s responsible for the 
synchrotron emission (leptonic jet models; for a recent review see, e.g., iBottcherl l2007al ) . or 
due to hadronic processes initiated by relativistic p rotons co-accelerated with the electrons 
(hadr onic models, for a recent discussion see, e.g., 
2003h . 



Miicke fc Protheroe 2001: Miicke et al. 



The quasar 3C279 {z = 0.536) is one of the best-observed flat spectrum radio quasars, 
in part because of its prominent 7-ray flare shortly after the launch of the Compton Gamma- 
Ray Observatory ( CGRO) in 1991. It was persistently detected by the Energetic Gamma-ray 
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Experiment Telescope (EGRET) on board CGRO each time it was observed, even in its very 
low quiescent states, e.g., in the winter of 1992 - 1993, and is known to vary in 7-ray 
flux by roughly two orders of magnitude (IMaraschi et al.l 1 1994 IWehrle et al.l Il998l ) . It has 
been monitored intensively at radio, optical, and more recently also X-ray frequencies, and 



has been the subjec t of intensive multi wavelength campaigns (e.g.. iMaraschi et al.l Il994 



Hartman et al.lll996l : IWehrle et al.lll998l ). Observations with the International Ultraviolet 
Explorer in the very low activity state of the source in December 1992 - January 1993 
revea led the existence of a thermal emission component, possibly related to an accretion 
disk JPian et al.lll999h . 



A complete compilation and modeling of all available SEDs s i multan eous with the 11 
EGRET observing epochs has been presented in iHartman et al.l (l2001a( ). The modeling 
was done using the time -dependent leptonic synchrotron self-C o mpton (SSC) -|- E x terna l 
Compton (EC) model of iBottcher. Mause. fc Schhckeiseil Jl997h : ISottcher fc BloonJ J2000h 
and yielded quite s atisfactory fits for all epochs. The results were consistent with other model 
fitting works (e.g., lBednareklll998l : ISikora et al.ll200ll : iModerski et al.l 120031 ) concluding that 
the X-ray - soft 7-ray portion of the SED might be dominated by SSC emission, while 
the EGRET emission might require an additional component, most likely external- Compton 
emission. 

During a recent observing carn paign by the Whole Earth Blazar Telescope (WEBT) 
collaboration (IBottcher et al.ll2007bl ) in the spring of 2006, intensive monitoring by the Ma- 
jor Atmospheric Gamma-Ray Imaging Cheren kov Telescope (MA GIC) yielded a positive 
detection at > 100 GeV on February 23, 2006 JAlbert et allboosh . This makes 3C279 the 
first quasar and (as of April 2009) the most distant object detected in VHE 7-rays. In this 
paper, we present the optical (BVRI) and X-ray (RXTE) data taken simultaneously with 
the MAGIC detection, and discuss the implications of this detection for current standard 
blazar jet models. 

Throughout this paper, we refer to a as the energy spectral index, Fj, [Jy] cx; z^^". A 
cosmology with = 0.3, Qa = 0.7, and Hq = 70 km s~^ Mpc~^ is used. In this cosmology 
the luminosity distance of 3C 279 at a redshift of z = 0.536 is = 3.08 Gpc. 



2. Observations and results 

3C 279 was observed in a WEBT campaign at radio, near-IR, and optical frequencies, 
throughout the spring of 2006. Details of the observations, data analysis, and implications 
of the optical variability patterns observed during that campaign have been published in 
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Fig. 1. — Light curves at radio (bottom), optical (center), and X-ray frequencies over the 
course of the WEBT campaign in the spring of 2006. The red points (right axis labels) 
indicate the energy spectral index in the 2-10 keV range. The vertical dashed line marks 
the day of the MAGIC VHE 7-ray detection. 
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Bottcher et al.l (l2007bl ). The source was simultaneously monitored with 3 pointings per 
week with the Rossi X-ray Timing Explorer (RXTE) Proportional Counter Array (PCA). 
We obtained the X-ray flux measurements with the PCA detector PC U2, using typical 



expos ure times of 2 ks for each pointing. The data reduction is described in IChatterjee et al. 
J2008h . 



Fig. [T] shows the radio, optical and X-ray light curves of 3C279 during spring 2006, 
along with the 2-10 keV energy spectral index as a function of time. The dashed vertical 
line marks the day of the MAGIC > 100 GeV 7-ray detection. While the source was overall 
in an extended optical high state {R ~ 14.5), no extraordinary variability in any optical 
(BVRI) band was observed at the time of the MAGIC detection. 

During most of December 2005 and January 2006, the X-ray flux of 3C279 was in a low 
state, near its historical minimum. Around Jan. 25, however, the source made a transition 
to a higher X-ray flux state with substantial variability in flux and spectral index on a 
characteristic time scale of ~ 10 days. The average flux increased to about a factor ~ 2 - 3 
compared to the low state. In the high state, there is a clear correlation between X-ray flux 
and spectral hardness, with the spectrum becoming harder as the flux increases. Statistical 
uncertainties preclude any conclusions about a fiux-hardness correlation in the low X-ray 
state. The VHE flare observed by MAGIC precedes an X-ray outburst wit h the highest X- 



ray fl ux measured since the major optical/X-ray outburst in 2001 (see, e.g. IChatterjee et al. 
2008h . by ~ 5 - 7 days. 



Fig. [2] compares several historical SEDs of 3C279 to the one measured on February 
23, 2006, along with the MAGIC VHE detection. The X-ray flux is comparable to the one 
observed during the major EGRET-detected 7-ray outburst in June 1991. The MAGIC data 
points show the measured flux, corrected for intergalactic 77 absorption, using the lowest 
plausible lev el of extragalactic ba ckground light ( EEL) cosmic i nfrare d light, according to 
the model of Primack et al. J2005h . As discussed in [Albert et al.l (120081 ) . this yields a best-fit 
energy index to the corrected VHE spectrum of avHE = l-9±0.9stat±0.5syst- Other currently 
discussed EBL models would predict a higher intri nsic VHE 7-ray flux and subst antially 
harder spectrum. In particular, a high-EBL model (IStecker. Malkan. &: Scullyll2006l ) would 



lead to an intrinsic VHE spectral index of avHE = —0.5 ± 1.2sta,t =1 = 0-5sv st. As we will 



discuss below, the VHE spectrum corrected with the iPrimack et al.l (120051 ) model already 
poses severe constraints on currently discussed blazar models. Those constraints would be 
even more restrictive and problematic for a higher EBL model. In particular, the high-EBL 
spectral index quoted above would locate the high-energy peak of the SED of 3C 279 at TeV 
energies, which would lead to even more unrealistic parameter choices for the leptonic models 
considered in the following section. We will therefore restrict our more detailed discussion 
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Fig. 2. — Compilation of broadband spectral energy distributions of 3C279, including the 
day of the MAGIC V HE detection, Febru ary 23, 2006. Data for PI and P2 are from 
f|2001al ). 2 003 data are fr o m ICollmar et al.l (j2004j ). and 2006 data are from 



Hartman et a. 



Collmar et al 



(120071 ) and lBottcher et al.l (!2007bl ). 
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to the VHE spectrum corrected by the low EBL modeL The optical spectrum, while clearly 
in an elevated state, shows about the same, steep spectral index Oopt ~ 1-7 as during lower 
optical flux states, indicating an underlying nonthermal electron spectral index of p = 4.4. 

The slopes of the radio and optical spectra indicate that the synchrotron peak was in 
the usual range where it has been observed in many previous observing campaigns, i.e., in 
the infrared regime, around i^sy ~ 5 x 10^^ Hz, corresponding to a dimensionless photon 
energy egy = /iz/sy/(meC^) ~ 4 x 10^''. This is consistent with the location of the sy nchrotron 



peak in a compilation of simultaneous multiwavelength data in mid-2006, shown by lMarscher 



( 120081 ). which included infrared coverage by the Spitzer Space Telescope. For the purpose of 
a quantitative analysis, the synchrotron peak flux may be estimated to be of the order of 
j/^sy ^ ]^gi3 jy Equally, the X-ray spectrum shows a quite typical shape as observed 
in previous high states of 3C 279, in particular the PI SED shown in Figs. [2]and|H This 
suggests that the X-ray - GeV 7-ray spectrum is similar to previously observed high states 
during the EGRET era. 

If a one-zone leptonic jet model (as discussed in the following section) applies, the VHE 
spectrum is expected to be at least as steep as the optical (synchrotron) flux. The spectral 
indices are expected to be similar if the 7-ray emission is produced by Compton scattering 
in the Thomson regime. If Klein-Nishina effects are important in the production of VHE 
7-rays, the resulting VHE spectrum would be even steeper than the synchrotron spectrum. 
As already indicated above, this would be in direct conflict with the observed relatively hard 
intrinsic VHE spectrum, even when corrected with a low EBL model. Therefore, in order not 
to predict a GeV 7-ray flux greatly in excess of any archival EGRET flux, it is reasonable to 
assume a 7-ray peak at ~ 10^^ - 10^^ Hz, corresponding to ~ 10^. The 7-ray peak flux is 
then uFJ ~ 5 x 10^^ Jy Hz. Previous modeling works of the SEDs of 3C 279 placed the peak 
of the high-energy component typically at frequencies around ~ 10^'^ Hz, i.e., about 1-2 
orders of magnitude lower than our new estimate, taking the MAGIC results into account. 
It is this shift of the inferred high-energy peak which will lead to quite dramatically different 
model implications for both leptonic and hadronic models compared to previous modeling 
efforts. 

From the spectral upturn in the UV in the P2-spectrum in Fig. [21 we can estimate a 
thermal (external) photon source with a luminosity of ~ 2 x 10^^ erg s~^, peaking at 
ud ~ lO^^ Hz (en ~ 10"^). 



- 8 - 



3. Implications for Leptonic Jet Models 

In this section, we consider wlietlier a one-zone leptonic blazar jet model can account 
for the February 23, 2006, SED of 3C279. We consider a scenario in which a nonthermal 
population of ultrarelativistic electrons produces, at the same time, the synchrotron emission 
from radio through UV and the 7-ray emission via Compton scattering of soft photons off the 
relativistic electrons. In general, we will assume that electrons are accelerated into a power- 
law distribution in electron energy, Q{'~f) = Qol^'^, in the range 71 < 7 < 72- The emission 
region has a radius Rb = lO^^-Rie cm. We define an escape time scale parameter ?7esc such 
that the escape time scale for relativistic electrons is Tosc = ^7osc Rb/c The interplay between 
radiative cooling and escape leads to the development of a spectral break in the electron 
spectrum at a Lorentz factor 7^, where the radiative cooling time scale equals the escape time 
scale. If the primarily injected electron distribution has a low-energy cutoff below the break 
Lorentz factor 7f, (the slow cooling regime), the spectral index of the electron distribution 
at 7 < 7fe is p = g, while above it the spectrum is steepened to p = q + 1. However, given 
the steep spectral index of the optical spectrum, implying p = 4.4 in the electron energy 
range synchrotron-radiating in the optical regime, a spectral break from p = 3.4 to p = 4.4 
would not produce a peak in the uF^ spectrum at the characteristic synchrotron frequency 
corresponding to 7?,. Therefore, it is more likely that the injected electron distribution has 
a high low-energy cutoff at 71 > In this case, particles at energies 7 < 71 result only 
from radiative cooling from higher energies, resulting in an electron spectrum with a spectral 
index p = 2 in the range 76 < 7 < 71, while above 71, we have p = g + 1, as in the slow 
cooling case. 



3.1. SSC model 



In the past, VHE emission has only been observed in high-frequency peaked BL Lac 
objects. In that case, a model interpreting the 7-ray emission as synchrotron self-Compton 
(SSC) emission has proven to be very succ essful, although r e cent o bservations of rapid vari- 
ability on ti r ne sca les of a few minutes flAharonian et al.l (120071 ) for PKS 2155-304, and 
Albert et al.l (120071) for Mrk 50 1) are posing serious challenges to this interpretation (e.g. 
Finke. Dermer fc Bottcherll2008l ). and the VHE emission of the intermediate BL Lac object 
W Comae, recently detected by VERITAS (lAcciari et al.ll2008l ) is more plausibly explained 
by Comptonization of an external radiation field. We point out that previous modeling 
efforts on 3C279, prior to the MAGIC detection, have c oncluded that an EC c omponent 
is strongly preferred to explain the GeV 7-ra y detection (IHartman et al.l l2001al : iBednarek 
19981 : ISikora et al.ll200ll : iModerski et al.l 120031 ). However, since the MAGIC points provide 
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a yet unexplored new constraint on blazar models for 3C279, it is worthwhile to revisit the 
SSC hypothesis in this paper. 

Given the synchrotron origin of the low-frequency peak at egy ~ 4 x 10~^ and the SSC 
origin of the 7-ray peak at ~ 10^, the Lorentz factor of electrons 7p radiating at the 
synchrotron and SSC peaks, can be estimated from 

7p= /ix^l.6xl0^ (1) 

V ^sy 

At the same time, the synchrotron peak frequency is given by 

= 4.2 X 10Sp5G^/(l + ^) Hz (2) 

where Bq is the (co-moving) magnetic field in Gauss, and D = 10 Di = (r[l — /Srcos^obs])"^ 
is the Doppler enhancement factor. This yields an estimate of the magnetic field and the 
Doppler factor of 

Bg Dir^ 7 xlQ-^ (3) 

This indicates that such a scenario would imply unrealistically low magnetic fields compared 
to standard values of ~ 1 G found from SED modeling of 3C279 in other states, as well 
as other blazar-type quasars. Even other known TeV blazars, whose SEDs can usually be 
well represented with SSC models, typically require magnetic fields of -B > 0.1 G, several 
orders of magnitude above the estimate found here for 3C279. We therefore conclude that a 
one-zone SSC model for the SED of 3C279 on February 23, 2006, including the VHE 7-ray 
emission is very problematic. 



3.2. External Compton 

The leptonic external-Compton scenario is based on the assumption that photons from 
an external, quasi-isotropic radiation field with dimensionless photon energy are Compton- 
upscattered to the observed 7-ray photon energies. As mentioned earher, if Klein-Nishina 
effects become important in the production of the VHE emission, the resulting VHE 7-ray 
spectrum would be even steeper than the observed synchrotron spectrum, in contradiction 
with the observed VHE spectrum, even when corrected for the lowest plausible level of EBL 
77 absorption. We therefore consider here only the possibility that the VHE 7-rays in a 
leptonic scenario are produced by Thomson scattering. Soft photons can be upscattered 
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effectively in the Thomson regime at most up to energies = l/e^. This indicates that a 
photon field with a characteristic photon energy of the accretion disk field at e/5 ~ 10~^ can 
effectively serve as the seed photon field for upscattering to the observed > 100 GeV 7-rays. 
We assume that a fraction tblr = 10~^r_i of the accretion disk radiation is reprocessed 
in the broad line region, which is located at an aver age distance -Rr t.r = 0.1i?BLR,-i pc 



from the central engine. HST near-UV spectroscopy (jPian et al.l l2005i) indicates that the 
total luminosity of the BLR in 3C279 is L^lr ~ tblr-^t) ~ 2 x 10^^ erg/s, motivating the 
above scaling in terms of r_i. In the co-moving frame, the external photons will thus have 
a characteristic energy of e'^ = Fe^. In addition to Eq. ([2]), we now have an independent 
estimate for 7p, namely 



We can use Eq. ([2]) to obtain an estimate for the magnetic field: 

Bg= 1.8 xlQ-^ TID^\ (5) 
The energy density of external photons in the co-moving frame can be expressed as 

Ue^t^-r-^ • 6 

47ri?^LR c 

From the 7-ray dominance, L^/Lgy ~ "u'ext/w's ~ 5, we can then obtain an independent 
estimate of the magnetic field of 

Bg = 1.0 t'_{'T,R^I^ _,. (7) 
Combining the magnetic field estimates and ([7]), we find 

i?BLR,-i = 57rl(' (8) 



which is in drastic contrast to the estimate of iPian et al.l (120051 ) of -Rblr ~ 3 x 10 ^ pc 



Considering the peak level of the synchrotron flux, we can use Eq. (8) of iBottcher et al. 



( I2OO3I ) to relate the magnetic field in the emission region with the equipartition fraction 
cb = u'b/u'^, i.e., the ratio of co-moving energy densities in the magnetic field and the 
nonrelativistic electron population: 
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o -lO'^D-^f ^27 f-w ej? \ 



4 /2 ^2 \ 1/7 

(9) 



Setting this equal to the magnetic-field estimate ([5]) yields 



9.7 X IQ-^Rl^Tl 



(10) 



Consequently, if we choose a conventional value of the Lorentz factor F ~ 15, we find an 
uncomfortably low magnetic field of -B ~ 0.03 G, corresponding to cb ~ 1.7 x 10^^i?fg, i.e., 
a far sub-equipartition magnetic field. Such a situation would make jet confinement very 
problematic, and is in contradiction with model results for 3C279 in other observing epochs 
and for other quasar-type blazars in general, where magnetic fields of typically i? ~ 1 - a 
few G are inferred, in approximate equipartition with the relativistic electron population. 

Alternatively, forcing the system to attain approximate equipartition would require us 
to assume an uncomfortably high Lorentz factor of F ~ 140 i?^g . This choice of a bulk 
Lorentz (and Doppler) factor would imply B ~ 0.25 G, and a low-energy cut-off of the 
injected electron population at 71 = 7p ~ 710 R^q . Apart from the fact that this is an order 
of magnitude larger than bulk Lorentz factors inferred from superluminal motion, it would 
require an implausibly close alignment of the jet with our line of sight, ^obs ~ 0.4°. We note 
that the magnetic-field estimate of Eq. ([5]) carries a proportionality B oc (e^/e^). Since our 
choice of ~ 10~^ is already close to the largest possible value to allow Thomson scattering 
to TeV 7-rays, the assumption of a different soft photon source (necessarily with a smaller 
eg) would worsen the problem of the unusually small inferred magnetic field. 

It has been noted by several authors that the 77 absorption of VHE 7-ray photons by the 
radiation field of the BLR may present another pro blem for a model of VHE 7-ray ernission 



inside the BLR of lumino us quasars in general fe.g. iDonea fc 



and 3C279 in particular (ILiu et al.ll2008l : ISitarek fc Bednarek 



rotheroell2003l : lReimerll2007r i 



20081 ). We therefore need to 



investigate the eff ects of 77 absorptio r i in t he BLR radiation field for the parameters we 
inferred above. In lBottcher &: Dermerl (jl995l ). the time- dependent 77 absorption signatures 
of an accretion disk fiare (refiected by BLR clouds) on VHE 7-ray emission were investigated. 
We have modified their approach for our purposes, adopting non-variable accretion- disk 
emission. The standard optically-thick accretion-disk spectrum is approximated by a spectral 
shape oc e^^'^e~^^^°, where Qd ~ 10~^ is the dimensionless inner disk temperature, 
©D = kTD i^/^meC^). Fig. [3] illustrates the dependence of the 77 absorption depth as a 
function of the dimensionless photon eriergy e and location of the 7-ray production site. Our 
results confirm the findings of iLiu et al.l fl2008n: With standard values of the BLR parameters 
, tblr ~ 0.1, -Rblr ~ 0.03 pc, as inferred by lPian et al.l (120051 ) . VHE 7-rays produced within 
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E/(m^ c ) 



Fig. 3. — Opacity for VHE 7-ray photons due to 77 absorption on the BLR radiation 
field. The labels denote the location of the 7-ray emitting region. Other parameters: = 



2 X 10^5 erg/s, 9^ = 10 



0.1 . Heavy (black) curves refer to the value of i?BLR,i 



0.03 pc as inferred by lPian et al.l (120051 ) with the outer edge of the BLR, -ReLR.out = 0.031 pc; 
light (blue) curves refer to -RBLR,in = 5.7 pc as inferred from Eq. [8] and -ReLR.out = 5.8 pc. 
The photon energy E is in the stationary AGN rest frame. 
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the BLR of 3C279 suffer severe 77 absorption by the same photon field that would serve as 
seed photon field for Compton scattering in a leptonic model. For the extreme parameters 
of r ~ 140, requiring -Rblr ~ 5.7 pc, tblr ~ 0.1, 77 absorption would hardly be a problem 
even out to multi-TeV 7-ray energies, if the 7-rays are produced close to the inner boundary 
of the BLR. Any level of 77 absorption in the BLR radiation field would imply an even 
higher level of intrinsic VHE 7-ray emission and therefore make the requirements for jet 
parameters in a one-zone leptonic model even more extreme. 

The red, solid curve in Fig. H] shows a leptonic model calculation with parameters similar 
to the quasi-equipartition case outhned above. We used a n equilibrium version of the time- 



dependent SSC -|- EC model of iBottcher fc Chiangl (120021 ). While the optical (synchrotron) 



spectrum and the level of the VHE 7-ray flux can be reproduced reasonably well, it is obvious 
that the X-ray flux is grossly underproduced. This is a consequence of the required, rather 
large low-energy cutoff at 71 ~ 700. The choice of a substantially lower 71 would extend 
both the synchrotron and the 7-ray spectra towards lower frequencies along the slopes of 
the optical and VHE 7-ray spectra and would therefore produce unreasonably large infrared 
and MeV - GeV 7-ray fluxes. In addition, it would require a larger jet power and therefore 
drive the system further out of equipartition (towards far sub-equipartition magnetic fields). 
Also, even though most of the > 100 GeV flux is produced by Compton scattering in the 
Thomson regime, the resulting VHE 7-ray spectrum appears to be steeper than the observed 



one, even with correction for the low EBL level according to iPrimack et al.l (120051 ). 



We therefore conclude that both the SSC and the external- Compton scenario for a one- 
zone, homogeneous jet model face severe problems representing the simultaneous SED of 
3C279 on February 23, 2006, including the VHE 7-ray emission. 

As noted earlier, the VHE 7-ray flare was not accompanied by any remarkable optical 
variability. This may be another hint that in 3C279 the optical and VHE 7-ray fluxes may 
be produced in separate emission regions. However, the calculation of the 77 opacity above 
indicates that even an inhomogeneous leptonic jet model would face severe problems in an 
external- Compton scenario. The dashed red curve in Fig. H] indicates that an SSC model 
can successfully reproduce the X-ray - VHE 7-ray spectrum with reasonable parameters, but 
fails to reproduce the optical spectrum. This may indicate support for a multi-zone leptonic 
model. This would require that the VH E 7-ray emission is produced far outside the BLR, 



possi bly in an internal shock scenario (jSpada et al.l 1200 ll : ISokolov. Marscher fc McHardy 



20041 ). This is indicated both by the 77 opacity argument as well as the low required 
magnetic field for the SSC fit presented in Fig. HI Such a model would imply that the X-ray 
emission is produced by low-energy electrons which have undergone substantial radiative 
cooling. One would therefore predict a delay of X-rays with respect to 7-rays on the order 
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v[Hz] 

Fig. 4. — Spectral fits to the SED of 3C279 on February 23, 2006: (solid [red]) using a leptonic 
External-Compton model with parameters similar to those derived in §3.2t (short-dashed 
[red]) leptonic SSC model fit only to the X-ray - 7-ray spectrum. Relevant parameters: 
7i = 10^, 72 = 10^ q = 2.3, Lj,e = 2.2 x 10^^ erg s-\ T^esc = 80, Rb = Q x 10^^ cm, 
r = D = 20, B = 0.2 G; (dot-dashed [maroon]) fit with the hadronic synchrotron-proton 
blazar model with internal (synchyrotron) photons only as targets for p7 pion production, 
and (long-dashed [maroon]) with synchrotron + external (BLR) photons as targets for p7 
pion production. See table [T] for parameters. 
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of the electron cooling time scale, 



obs 
'cooLSSC 



rrtpC 



D{A/?,)caTu'B{Lssc/L,y)l ^^^^ 
Using the value of 5 = 0.2 G used for the model shown in Fig. HI and scaling the electron 
energy 7 = 10^73, we find the expected 7-ray vs. X-ray delay as r ~ 277^^ hr. We note, 
however, that an intensive search for inter-band ti me delays betw e en opt ical, X-ray and 7- 
ray emission of 3C279 during the EGRET era by iHartman et al.l (j2001bl ) did not find any 
evidence for systematic 7-ray vs. X-ray delays. 



Hadronic Models 



If relativistic protons at energies above the threshold for p7 pion production are present 
in the jet of 3C 279, hadronic interactions must be considered in blazar emission models. 
For the present modeling w 'e use the hadronic Synchrotron-Proton Blazar (SPB) model of 
Miicke fc Protherod (120011 ): iMiicke et al.l (120031 ). In its original version, this model was best 
suited for X-ray selected BL Lac objects with very weak or absent external radiation fields. 
However, the quasar 3C 279 is known to have strong accretion disk and BLR line emission 
(see §3.21) . We therefore extend the (one-zone) SPB model to account also for target photon 
fields external to the jet. 

Relativistic electrons (e) and protons (p) with a power law index qe = qp are injected 
instantaneously into a spherical emission region, or blob, which is moving with relativistic 
velocity along the jet axis. In the strongly magnetized (with field strength B = const.) blob, 
the primary electrons lose their energy predominantly through the synchrotron channel. The 
resulting synchrotron radiation generally dominates the low energy component of the blazar 
SED, and serves as target photon field for proton-photon interactions and pair (synchrotron) 
cascading. In our extension of the SPB model we add th e photon field from the BLR, 
simp lified as an isotropic distribution in the jet frame (see iReimer. Bottcher fc Postnikov 



2005 



for a discussion of this approximation) as an additional target for particle-photon 
interactions and cascading. We properly take into account the Doppler boost of the average 
photon energy and photon energy dens ity of the external photon fiel d. In the present work 
we use the approximation advocated by lTavecchio &: Ghisellinil ( 120081 ) for the BLR radiation 
field: In the co-moving jet frame a blackbody spectrum with temperature T' cx Tu^ya {vhya 
the Lya line frequency) is used to approximate the bo osted BLR em ission from the stationary 
AGN to the blob frame. Local absorption of 7-rays (IReimerl 120071 ) external to the emission 



region is taken into account as well. The injected relativistic protons suffer energy losses from 
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photomeson production, Bethe-Heitler pair production, synchrotron radiation and adiabatic 
expansion. Charged particles produced as secondaries in the photomeson production channel 
may suffer synchrotron losses as well prior to their decay. This is particularly relevant for 
charged pions and muons (henceforth named /x-synchrotron radiation). All high energy 
photons may initiate pair cascades, which redistributes a fraction of the photon power from 
high to lower energies where the photons can eventually escape the emission region. 

In the framework of the hadronic SPB model the injection electron spectrum is primarily 
constrained by the optical and radio data. For = 2.1- 2.2 the primary electron synchrotron 
spectrum above the synchrotron self-absorption turnover shows a rather flat spectrum from 
the injected particle spectrum, modified by synchrotron losses, followed by a steep tail due 
to the cutoff of the electron distribution at particle Lorentz factors ~ 10'^. The high energy 
component of the SED of 3C 279 is constrained by the RXTE and MAGIC data. In the SPB 
model the RXTE data can in general be explained by either proton synchrotron radiation, or 
a reprocessed/cascade component. The former suggests, however, extremely long loss time 
scales (of order years for typical field strengths and Doppler factors in the SPB model), which 
is difficult to reconcile with the observed day-scale variability. We therefore concentrate on 
the second option of reprocessed radiation dominating the X-ray band. This picture is also 
strengthened by the softness of the X-ray spectrum, which may indicate the appearance of 
reprocessed (through 77 pair production) radiation in this energy range. 

The left panel in Fig. O shows a typical model fit to the SED of 3C 279 of February 
23, 2006, for the case of a negligible external photon field at the location of the 7-ray 
emission region. While radiation in the 7-ray band is dominated by proton and n synchrotron 
emission, reprocessed proton synchrotron radiation determines the photon emission in the 
X-ray band. For the required bulk Doppler factors D = 10 - 14, and assuming a size of the 
emission region of order ~ lO'^^ cm, the energy density of the internal jet target photon field 
amounts to u'sy ~ a few 10^° - 10^^ eV cm~^. With these values, a delay between the TeV 
and X-ray band of a few days may be explainable. The strong magnetic field strengths of 
40 - 60 G imply losses due to proton synchrotron radiation in the TeV band on hour time 
scales. Interestingly, all models representing the simultaneous data of February 23, 2006, 
reasonably well require a cutoff of the injected proton spectrum at a few 10^ GeV, which 
is significantly lower than what is needed for HBL-like TeV-blazars. The injected proton 
energy density u'p is of the order 10^ — 10^ erg cm~^, not too far from the equipartition value 
^^p gq^j (here u'p ^ 2 - 17^^p gq^j), and a total jet luminosity (as measured in the galaxy rest 
frame) of order io^^~^^ erg s~^. 

If the gamma-ray emission region is rather close to the BLR, the external photon field 
has to be taken into account as an additional target for particle-photon interactions and pair 
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Table 1. Parameters for the hadronic synchrotron-proton blazar model fits to 3C279 
(Figs. HlandE]) for (2nd col.) internal synchrotron photons only, and (3rd col.) internal 
synchrotron + external (BLR) photons as targets for p7 pion production. 



Parameter 


symbol 


SPB (sy.) 


SPB (sy. + ext.) 


Jet power 


Lj [erg 


1 X 10*9 


6.6 xlO**** 


Doppler factor 


D 


12 


12 


Blob radius 


Rb 


3 X IQi'^ cm 


3 X IQi" cm 


Maximum proton energy 




5 xlO^ GeV 


4.5 X 10^ GeV 


Magnetic field 


B 


40 G 


60 G 


Particle spectral index 


qp = Qe 


2.2 


2.2 


Proton energy density 


Up [erg cm~'^] 


1.2 x 10^ 


540 


Electron/proton density ratio 


rie/rip 


2 X 10-3 


3 X 10-3 


Temperature of BLR emission 






3 eV 


BLR photon energy density 


mblr 




10~3 ergs cm-3 


BLR radius 


^in.BLR 




0.1 pc 


Location of emission region 


Rem 




0.18 pc 




Loq (v/Hz) Loq (iz/Hz) 

Fig. 5. — Hadronic fits to the SED of 3C279. Left: Pure synchrotron-proton blazar model 
with only intrinsic synchrotron photons as targets for p7 pion production; right: synchrotron 
proton blazar model including external photon field from the BLR. See table[I]for parameters. 
Individual radiation components are: dashed: proton synchrotron and cascade; triple-dot- 
dashed: synchrotron and cascade; dotted: 7r° cascade; dot-dashed: n"^ cascade. 
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cascading. The observed line and disk luminosity of 3C 279 implies tblr = Lblr/Ld ~ 0.08. 
The additional target photon field in the optical/UV band leads to enhanced reprocessing 
in the hadronic model, further softening the spectrum in the X-ray band. Signatures of 7- 
ray absorption due to the additional narrow-banded external photon field are clearly visible 
in the MAGIC energy range. The right panel in Fig. [5] shows an example of a model fit 
representing again the February 23, 2006, data. Proton synchrotron radiation dominates in 
the 7-ray band. The reprocessed component is more pronounced here than for the case of 
internal target photon fields only, which yields a more appropriate description of a potentially 
steep X-ray spectrum. 77 absorption of VHE 7-ray photons in the external radiation field 
inevitably leads to a steepening of the VHE 7-ray spectrum. The fit shown in Fig. [Sjo 
presents an acceptable representation of the MAGIC spectrum. Should a higher EBL level 
imply a harder intrinsic VHE 7-ray spectrum, this fit could be modified by choosing a lower 
energy density of the external radiation. Parameter values for hadronic model fits are given 
in Table [H For the fit with external photons, the injected proton energy densities are at 
most a factor 4 above the equipartition value and the required total jet luminosity of ~ (1 
- 7) X 10^^ erg s~^ is somewhat lower than for the case of internal target photon fields only. 



5. Summary 

We have presented simultaneous optical and X-ray spectral information to the recent 
MAGIC detection of the quasar-type blazar 3C279 in February 2006. The source was shown 
to be in an elevated optical state, but showed no substantial optical variability and a rather 
steep optical spectrum during the MAGIC detection. The GeV - TeV fiare preceded an X-ray 
fiare by about 5-7 days. We have presented the simultaneous broadband (optical - X-ray 
- VHE 7-rays) SED of 3C279 and discuss its implication for one-zone leptonic jet models. 
We found that an SSC model is extremely problematic as it would imply unreasonably low 
magnetic fields. Also an external-Compton interpretation has problems with an unusually 
low magnetic field of -B ~ 0.03 G, implying an equipartition ratio of ~ 10~^. Alternatively, 
approximate equipartition can be achieved with a bulk Lorentz factor of F ~ 140 -R^e^''^, which 
appears equally unlikely. These constraints have been inferred for a correction of the VHE 
spectrum for 77 absorption by the lowest plausible EBL model and without taking into 
account internal 77 absorption by photons from the accretion disk or the broad-line region. 
Any higher EBL level or a substantial amount of internal 77 absorption would lead to even 
more extreme constraints. 

We therefore conclude that a simple homogeneous, one-zone leptonic jet model has 
serious problems reproducing the SED of 3C279 on February 23, 2006, which includes the 
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recent VHE 7-ray detection by MAGIC. The lack of correlated optical - 7-ray variability 
suggests, instead, a multi-zone model in which the optical and 7-ray fluxes are produced in 
separate regions along the jet. We have shown that a leptonic SSC fit to the X-ray — VHE 
7-ray spectrum alone can be achieved with parameters quite typical for quasar-type blazars. 

Alternatively, the hadronic synchrotron-proton blazar model is able to provide an ac- 
ceptable fit to the SED of 3C279. Both a pure synchrotron-proton blazar model (without 
external photons as targets for p7 pion production) and a model with a substantial con- 
tribution from external target photons can reproduce the observed SED up to VHE 7-ray 
energies very well. If no external photon field is included in the model, the relevant proton 
synchrotron energy loss time scale is of the order of years and would therefore be inconsis- 
tent with the observed day-scale variability if proton synchrotron radiation was dominant 
in the X-ray regime. This model therefore requires a high internal radiation energy density 
in order for proton-synchrotron induced cascades to dominate the X-ray emission. Alter- 
natively, an additional target photon contribution from external sources (in particular, the 
BLR) is able to overcome this problem, and such a model seems to provide an appropri- 
ate fit to the observed SED of 3C279 from optical to VHE 7-rays. However, both versions 
of the hadronic synchrotron proton blazar model require a rather extreme jet power of 
Lj ~ 10^^ erg s^^. For comparison, a multi-zone leptonic model requires a jet power in 
leptons alone of Lj^e ~ 2.2 x 10^^ erg s~^. When including an equal number of cold protons, 
the total jet power requirement increases to -£/j,tot ~ 2.6 x 10^^ erg s~^. 

We point out that our conclusions arc based on rather incomplete frequency coverage 
during the MAGIC detection, and rely on the apparent similarity of the SED, at least from 
radio through X-rays, with previously observed high states of 3C279 during the EGRET era. 
Future simultaneous Fermi LAT and TeV 7-ray observations (with a detection at > 100 GeV 
7-rays) would provide crucial tests to our conjectures put forth in this paper. 

The work of M. Bottcher and was partially supported by NASA through XMM-Newton 
GO grant awards NNX07AR88G and NNX08AD67G. The work of A. Marscher was funded in 
part by the National Science Foundation through grant AST-0406865 and by NASA through 
RXTE Guest Investigator grant NNX06AG86G and Astrophysical Data Analysis Program 
grant NNX08AJ64G. 



REFERENCES 

Acciari, V. A., et al, 2008, ApJ, 684, L73 



- 20 - 

Aharonian, F. A., et al., 2007, ApJ, 664, L71 

Albert, J., et al., 2007, ApJ, 669, 862 

Albert, J., et al., 2008, Science, vol. 320, no. 5884, p. 1752 

Bednarek, W., 1998, A&A, 336, 123 

Bdttcher, M., & Dermer, C. D., 1995, A&A, 302, 37 

Bottcher, M., Mause, H., & Schlickeiser, R., 1997, A&A, 324, 395 

B6ttcher, M., & Bloom, S. D., 2000, AJ, 119, 469 

Bottcher, M., & Chiang, J., 2002, ApJ, 581, 127 

B6ttcher, M., et al., 2003, ApJ, 596, 847 

Bottcher, M., 2007a, in proc. "The Multimessenger Approach to Gamma-Ray Sources", 
ApSS, 309, 95 

Bottcher, M., et al., 2007, ApJ, 670, 970 

Chatterjee, R., et al., 2008, ApJ, 689, in press (larXiv:0808.2194p 



CoUmar, W., et al., 2004, in proc. "The INTEGRAL Universe", Eds. G. Lichti, V. 
Schonfelder, & C. Winkler, ESA-SP 552, p. 555 

Collmar, W., et al., 2007, in proc. "The Obscured Universe", Eds. S. Grebenev, R. Sunyaev, 
& C. Winkler, ESA-SP 662, p. 207 

Donea, A.-C., & Protheroe, R. J., 2003, Astrop. Phys., vol. 18, Issue 4, p. 377 

Finke, J. D., Dermer, C. D., & Bottcher, M., 2008, ApJ, 686, 181 

Hartman, R. C., et al., 1996, ApJ, 461, 698 

Hartman, R. C., et al., 1999, ApJS, 123, 79 

Hartman, R. C., et al., 2001a, ApJ, 553, 683 

Hartman, R. C., et al., 2001b, ApJ, 558, 583 

Liu, H. T., Bai, J. M., & Ma, L., 2008, ApJ, in press farXiv:0807.3133) 



- 21 - 



Marscher, A. P., in proc. of "Workshop on Blazar Variability across the Electromagnetic 
Spectrum" , Palaiseau, France, 2008, in press 

Maraschi, L., et al., 1994, ApJ, 435, L91 

Moderski, R., Sikora, M., Blazejowski, M., 2003, A&A, 406, 855 
Miicke, A., & Protheroe, R. J., 2001, Astropart. Phys., 15, 121 

Miicke, A., Protheroe, R. J., Engel, R., Rachen, J. P., & Stanev, T., 2003, Astropart. Phys., 
18, 593 

Pian, E., et al., 1999, ApJ, 521, 112 

Pian, E., Falomo, R., & Treves, A., 2005, MNRAS, 361, 919 

Primack, J. R., Bullock, J. S., & Somerville, R. S., in "High- Energy Gamma- Ray Astron- 
omy", AIP Conf. Series, F. Aharonian, H. Voelk, D. Horns, Eds. (AIP, Heidelberg, 
2005), vol. 745, p. 23 

Reimer, A., Bottcher, M., & Postnikov, S., 2005, ApJ, 630, 186 

Reimer, A., 2007, ApJ, 665, 1023 

Sikora, M., Blazejowski, M., Begelman, M. C, & Moderski, R., 2001, ApJ, 554, 1; Erratum: 
ApJ, 561, 1154 (2001) 

Sitarek, J., & Bednarek, W., 2008, MNRAS, in press (larXiv: 0807.4228]) 

Sokolov, A., Marscher, A. P., & McHardy, I. M., 2004, ApJ, 613, 725 

Spada, M., Ghisellini, G., Lazzati, D., & Celotti, A., 2001, MNRAS, 325, 1559 

Stecker, F. W., Malkan, M. A., & Scully, S. T., 2006, ApJ, 648, 774 

Tavecchio, F., & Ghisellini, G., 2008, MNRAS, 385, L98 

Wehrle, A. E., et al., 1998, ApJ, 497, 178 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



